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Abstract. We develop a physical model of the zodiacal cloud incor- 
porating the real dust sources of asteroidal, cometary, and kuiperoidal 
origin. Using the inferred distribution of the zodiacal dust, we compute 
its thermal emission and scattering at several wavelengths (1.25, 5, and 
20 pm) as a function of NGST location assumed to be at 1 AU or 3 AU. 
Areas on the sky with a minimum of zodiacal light are determined. 


1. Physical model for zodiacal thermal radiation and scattering 

Improvements in the zodiacal light emission and scattering to be given by space 
observations at 3 AU, compared with observations near the Earth, were dis- 
cussed by Mather and Beichman (1996). Unfortunately, a rather accurate multi- 
parametric model of the zodiacal brightness derived by Kelsall et al (1998) from 
the COBE data cannot be reliably extrapolated to heliocentric distances as large 
as 3 AU. We have developed a physical model of the zodiacal cloud incorporating 
the real dust sources of asteroidal and cometary origin, which makes it possible 
to evaluate quantitatively the zodiacal light emission and scattering throughout 
the Solar system ( Gorkavyi et al. 1997a). This model considerably improves 
our previous ‘reference model’ based on the use of the continuity equation for 
distribution function of dust particles ( Gorkavyi et al. 1997b ,c, 1998) and en- 
ables us to obtain more reliable results than a phenomenological modelling of 
the zodiacal light (e.g. Ebbets 1998). Below, we describe an improved model 
that represents a 3D-grid containing 45 x 180 X 244 = 2 * 10 6 cells with a step 
in (heliocentric latitude p , longitude A, and radius R) to be (2°, 2°, 0.025 R 
[AU]). Using a new numerical approach to the dynamics of minor bodies and 
dust particles, we increase the number of particle positions employed in each 
model to 10 11 without using a supercomputer ( Ozernoy et al. 2000 } Gorkavyi 
et al. 2000). We compute here the distribution of the zodiacal dust emission 
and scattering in the Solar system. The processes influencing the dust particle 
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dynamics include gravitational scattering on all planets, except Pluto (inclina- 
tions and precession of the planets are neglected), mean motion resonances, and 
the Poynting- Robertson/solar wind drags. 

We employ here 931 sources of dust particles (with eccentricities e < 1), 
which include 284 asteroids, 451 short-period comets, and 196 Kuiper belt ob- 
jects. We adopt a two-component approximation for dust particle size distribu- 
tion, which is characterized by the parameter p = 0.285 and 0.057 for (l-2)/im 
and (5-10)/xm particles, respectively. 

The inferred dust density distribution, n(R,<p , A), of the IPD at heliocentric 
distances 0.5 < R < 100 AU enables us to compute the scattered light and the 
thermal emission of the zodiacal cloud as a function of the observer’s latitude 
and longitude. 

The thermal emission and scattering of the zodiacal cloud are given by: 

/„ oc f n(R , <p t A) [(1 - A)B I/ (T) + AF ? *(©)] ds 

where B v is the Planck function at frequency v, A is albedo, T is the dust 
temperature, T(R.) = TqR~ ^ with To = 286°K and S = 0.467 ( Kelsall et al. 
1998), F® oc 1 /R 2 is the solar flux, $(0) is the phase function at scattering 
angle 0 ( Kelsall et al. 1998). 

Note that although the phase function only weakly depends on v, the value 
of I„ strongly depends on v through the solar flux. Besides, the Planck function 
strongly depends on heliocentric distance, especially in the near IR (1-5 \i m). 


2. The zodiacal light at R=1 and 3 AU 

Using the inferred distribution of the zodiacal dust, we have computed a variety 
of zodiacal light maps, both for thermal emission and scattered components, 
at different locations (R,Z) of the NGST assumed to be at 1 AU or 3 AU 
(see Fig. 1). The brightness as a function of latitude <p and longitude A (in 
the telescope’s frame) is given in logarithmic scale (the neighboring contour 
intensities differ by \/e). The Sun’s position is (0,0) for Z = 0 and is shifted to 
a negative <p for Z = 0.25 AU. 

At each location, there is a minimum in the zodiacal light which can be seen 
as a ‘dark spot’ (or several spots) in the computed Figures. The positions of 
those minima are explained by an interplay between dependencies of emissivity 
upon density and temperature. 

Our numerical results concerning the computed zodiacal light at different 
locations of NGST are summarized in Tables 1 and 2. The data presented in 
Table 2 for R = 1 AU, 2nd column can be compared with the COBE data 
/(0°,90°)//(90°) = 3.5 (1.25//m band), = 2.3 (5/*m band), and = 3.3 (20/xm 
band), which implies that the major contribution into the zodiacal light at 1 
AU is given by small asteroidal (al) and cometary (cl) dust components in 
approximately equal proportions. 
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Table 1. Components of the model zodiacal cloud. 


Origin/Size 
of dust particles 

# of sources® 

# of particles 
with e < t at t = 0 ° 

# of particles 
near Earth® 

# of positions 
in our 3D-modei 

asteroidal/l^m (al) 

110+110 

110+110 

57+85 

7 • 10 9 

cometary/X/im (cl) 

2128+1X2 

276+112 

73+78 

7- 10 9 

kuiperoidal/l^xm (kl) 

100+100 

96+100 

11 + 12 

6 • 10 10 

asteroidal/5/xm (a5) 

32+32 

32+32 

32+32 

2 • 10 10 

cometary/5^m (c5) 

61+3 

60+3 

40+2 

1 • L0 l ° 


° Pericentral start + apocentral start 


Table 2. Brightness of different components of the zodiacal cloud 

assuming the intensity of each component near the Earth be I(y> = 90°, A) = 1 at the pole 


0* 

R = 1 AU, 2 

2 b 3 c 

' = 0 AU 
4 d 5* 

R = 1 AU, Z 
1» 2 b 3 C 

= 0.25 AU 
4 d 5 e 

1° 

R = 
2 b 

3 AU, Z = 0 AU 
3 C 4 d 5' 







1.25 /im band 







al 

5.0 

3.4 

.93 

75 

.18 

1.7 

1.3 

.17 

87 

.068 

.30 

.12 

.064 

74 

cl 

2.1 

1.4 

.92 

75 

.38 

2.0 

1.5 

.38 

89 

.057 

.14 

.092 

.054 

54 

kl 

4.1 

2.8 

.95 

73 

.16 

2.0 

1.8 

.15 

89 

.049 

.57 

.62 

.046 

87 

a5 

7.5 

5.1 

.90 

75 

.064 

1.1 

.93 

.063 

88 

.018 

.060 

.019 

.012 

45 

c5 

1.8 

1.2 

.91 

75 

.49 

1.6 

1.1 

.49 

89 

.045 

.077 

.052 

.038 

57 








5 fj . m band-f 







al 

3.0 

1.0 

.80 

41 

.17 

.79 

.31 

.16 

82 

.70 

1.9 

.50 

.47 

38 

cl 

1.4 

.46 

.46 

0 

.36 

1.2 

.42 

.31 

65 

.53 

.69 

.21 

.21 

0 

kl 

2.4 

.78 

.75 

39 

.15 

.79 

.29 

.14 

89 

.44 

.73 

.22 

.20 

10 

a5 

5.0 

1.8 

.87 

67 

.062 

.51 

.22 

.058 

81 

.19 

.55 

.12 

.11 

45 

c5 

1.3 

.43 

.43 

14 

.47 

1.1 

.37 

.30 

14 

.42 

.51 

.14 

.14 

0 







20^m band 







al 

5.3 

2.8 

.98 

74 

.20 

1.9 

1.2 

.20 

82 

.10 

.45 

.18 

.098 

74 

cl 

2.1 

1.1 

.95 

73 

.43 

2.1 

1.1 

.43 

89 

.083 

.18 

.089 

.071 

46 

kl 

3.9 

1.9 

.94 

83 

.17 

1.8 

1.1 

.16 

89 

.068 

.36 

.24 

.063 

78 

a5 

8.8 

4.5 

.99 

83 

.080 

1.3 

.82 

.079 

88 

.029 

.096 

.025 

.019 

45 

c5 

1.8 

.87 

.83 

16 

.56 

1.7 

.81 

.56 

89 

.066 

.10 

.042 

.042 

0 


‘Type of dust listed in Table 1. 
intensity at the north pole, I po u 

6 intensity at the ecliptic plane in the direction of 90° from the Sun, /(y? = 0°, A = 90°) 
intensity at the ecliptic plane in anti-sun direction, /(y? = 0°, A = 180°) 
d / min (y?min, A = 180°), where y? ra in is the latitude at which the zodiacal light brightness 
is minimal in the anti-sun direction (A = 180°) 

a y>min [°] 

f All 5^m-intensities at R = 3 AU, Z = 0 need to be multiplied by 10“ 3 


3. Conclusions 

1. The structure of the zodiacal dust cloud computed in the present work is 
substantially non-uniform: 
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(i) near the Earth, the thickness of the dust layer is the largest for cometary 
particles and the smallest for asteroidal particles, with kuiperoidal particles oc- 
cupying an intermediate position; 

(ii) the larger the size of asteroidal dust particles, the thinner is the layer 
comprised of such particles; for cometary particles a reverse, although a weaker, 
dependence takes place; 

(iii) the partial contribution of particles of different origin and size changes 
with heliocentric distance: at R = 3 AU the contribution of asteroidal particles 
(especially of large size) into the zodiacal light emission is sharply decreased com- 
pared to the cometary component. At R = 3 AU, the contribution of kuiperoidal 
dust is substantial in the ecliptic plane, especially for the scattered light; 

(iv) The latitudinal dependence of the zodiacal emission is different for 
different components of the dust cloud. In the anti-sun direction (A = 180°), 
the maximum of the asteroidal dust emission at 5 fim is reached at the ecliptic 
plane, whereas the brightness of the cometary component has its minimum here. 

2. Observations made from R = 1 AU, Z — 0.25 AU at all wavelengths 
would detect the brightness of the zodiacal light in the direction of the pole at 
level of 25-30% of that seen near the Earth. 

3. At R = 3 AU, the zodiacal light brightness in the direction of the pole 
decreases to 5 - 6% at 1.25 /im; to 0.05 - 0.07% at 5 /im, and to 8 - 10% at 
20 pm compared to that seen near the Earth. 

4. At each wavelength, there are certain regions on the sky having a minimal 
zodiacal light brightness. At short wavelengths (~ 1.25 /xm), these regions are 
around the poles, and at long wavelengths (such as 5 to 20 /im) these regions 
are shifted toward the ecliptic. This result needs to be taken into consideration 
while planning observations with NGST. 

5. Further improvements in our physical modelling of the zodiacal light are 
feasible. Using COBE and IRAS data, the contribution of dust components of 
different origin into the zodiacal cloud could be determined. Those improve- 
ments would make far going implications possible for extragalactic astronomy 
and cosmology. 
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Figure 1. Emission of a small (r = 1 — 2/nn) asteroidal dust in the 
5/nrt band 
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